Topological constraints on magnetostatic traps 
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We theoretically investigate properties of magnetostatic traps for cold atoms that are subject to 
externally applied uniform fields. We show that loffe Pritchard traps and other stationary points of 
B are confined to a two-dimensional curved surface, or manifold M, defined by det{dBi/dxj) = 0. 
We describe how stationary points can be moved over the manifold by applying external uniform 
fields. The manifold also plays an important role in the behavior of points of zero field. Field zeroes 
occur in two distinct types, in separate regions of space divided by the manifold. Pairs of zeroes 
of opposite type can be created or annihilated on the manifold. Finally, we give examples of the 
manifold for cases of practical interest. 

PACS numbers: 03.75.Be, 41.20.Gz, 32.80.Pj 
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I. INTRODUCTION 

Magnetic trapping of neutral particles has first been 
achieved for cold neutrons [l| and has since become a 
widely used tool in cold-atom physics . More recently, 
the flexibility to design complicated magnetic trapping 
potentials has been boosted tremendously by the devel- 
opment of atom chips [^ H, p. The field sources are 
defined by microfabrication on a planar substrate, tak- 
ing the form of either current-carrying wires or patterns 
in a permanent magnetic film |y, III 131 • 

Magnetostatic traps are defined by magnetic field min- 
ima. In this paper we investigate the occurrence of field 
minima from a general perspective Q. We introduce a 
novel conceptual tool, a curved surface, or manifold A4, 
to which all stationary points of B (non-zero minima and 
saddle points) are confined. We derive expressions for the 
movement of stationary points over this manifold, in re- 
sponse to a change of an external uniform control field. 
We also show that the same manifold plays an important 
role in the creation and merging of field zeroes. 

The typical application that we have in mind is a situ- 
ation where a magnetic field configuration is fixed by e.g. 
permanent magnets and control of the field is limited to 
the application of uniform external fields. This situation 
occurs for instance in atom chip experiments [a, B 13; 
where field gradients can become very large. Control 
of the movement of loffe-Pritchard (IP) traps [lOl lll| is 
of importance in loading procedures and in experiments 
that require dynamical splitting or movement of atomic 
clouds. During loading for instance it is important to 
avoid regions of zero field, since this will lead to losses 
due to Majorana spin flips to untrappcd states. It is also 
important to avoid unwanted splitting of the trap during 
the transport and compression of the cloud to the final 
trap. Furthermore, quantum information processing ap- 



plications on an atom chip may require the movement of 
qubits to regions where they can be 'read out' or manip- 
ulated. In this case it is also of importance to keep track 
of the individual phase evolutions of the atoms, i.e. to 
control the trapping parameters during transport. 

Although we investigate here magnetic traps, most of 
our conclusions also apply to traps based on electrostatic 
fields insofar as they rely on the field being rotation and 
divergence free. Electrostatic traps can be used to trap 
molecules with an electric dipole moment [ij, [iJ, [3 • 

This paper is structured as follows. After introducing 
our notation in Sec. ^ in Sec. Illll we derive the expres- 
sion for the manifold to which stationary points must be 
confined. We also show how to create a IP trap in a 
given point on this manifold. In Sec. IIVI we derive an 
expression for the movement of stationary points along 
the manifold, under the influence of an external uniform 
field. In Sec. E] we investigate the relationship between 
the manifold and points of zero field. We show how field 
zeroes can be moved and how pairs of zeroes can be cre- 
ated and annihilated on the manifold. Finally, in Sec. I VII 
we investigate the shape of the manifold for some cases 
of experimental interest. 



II. NOTATION 

Magnetic traps are usually operated in the regime 
where moving particles experience a magnetic field that 
varies slowly compared to the Larmor spin precession fre- 
quency. The spin component parallel to the field is then 
conserved due to adiabatic following of the local direction 
of the magnetic field. The effective potential is propor- 
tional to the modulus of the field: B{r) = |B(r)|. We are 
interested in stationary points and trapping frequencies 
in this potential. For this purpose it is equivalent to use 
i3^(r), since a minimum or saddle point of B is also a 
minimum or saddle point of B^. We shall use B^{r) for 
convenience and define 
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U{r) ^B^{r). 



(1) 



Throughout this paper we adopt the convention that 
summation over repeated indices is impHed, e.g. BiBi = 
Bl + By + B^ ~ B"^. Points where U is stationary are 
defined by diU = dU/dxi = for all i = 1,2, 3. In order 
to decide whether a stationary point is a local minimum 
or a saddle point, wc will need also the second deriva- 
tives. Therefore, let us expand B(r) to second order in 
the relative coordinates Xi around some point of interest, 



B, 



gtjXj 



(-> CijkXjXk 



h.o.t. 



(2) 



where h.o.t. denotes higher order terms, gij = djBi is 
a tensor describing the gradient of the vector field B(r) 
and Cijk — Cikj = djdkBi is a curvature tensor. 

Using Maxwell's equations for stationary fields in vac- 
uum we can impose some restrictions on the tensor com- 
ponents gij and Cijk- From the conditions divB = 
and curlB = for stationary fields in empty space we 
see that the gradient tensor must be both traceless and 
symmetric. 



U, not of B. For stationary points at nonzero field, such 
as the minimum in a loffe-Pritchard trap, we must require 
that gip has an eigenvalue zero, i.e. 



det (gij) = 0. 



(8) 



Furthermore for a stationary point to occur, the field 
Ui must be parallel to the eigenvector of gip with eigen- 
value zero. In the case of a IP trap this direction is 
usually called the axial direction. The above condition, 
Eq. © , expresses the fact that a IP trap requires a point 
where the magnetic field locally looks like a cylindrical 
quadrupole field, and that the axis of the quadrupole is 
the trap axis. 

Since the gradient is a function of the spatial coor- 
dinates, gij = gij{v), the condition of Eq. (jSJ defines 
a two-dimensional curved surface which we shall call the 
manifold M . The points on the manifold are those points 
in space where a stationary point for B can be created by 
choosing Ui along the zero eigenvector of gij . Some exam- 
ples of such manifolds in situations of practical interest 
are shown in Sec. IVII 



9t3 





9ji 



(3) 

(4) 



B. Absence of field maxima in empty space 



This leaves five independent parameters for gij, which 
can be interpreted as follows. In the coordinate frame 
where gij is diagonal, two independent gradients can be 
chosen. Three angles are needed to specify the orienta- 
tion of the coordinate frame. 

Similarly, the curvature tensor Cijk must be fully sym- 
metric under permutation of indices and all its partial 
traces must vanish: 



^tij 



^iji 



^jii 







C-ijk — C-iUj — Ckji 



(5) 
(6) 



This leaves seven independent parameters for c^fc. 

Throughout the paper we adopt the convention that 
eigenvalues of a tensor are written in capital letters and 
eigenvectors are written in capital bold face letters. Thus 
(Gi, G2, G3) are the eigenvectors oigij and (Gi, 62, G3) 
are the corresponding eigenvalues. 



III. STATIONARY POINTS 

A. The manifold det{gij) = 

In order to find the stationary points of L'^ = BiBi, we 
substitute the expansion of Eq. ^ and collect terms up 
to second order in the coordinates. 



If the condition Uigip = is fulfilled, Eq. Q is simpli- 
fied (up to second order) to: 



U Uilli -i" TjkXjXk, 



where we defined 



f^jk — ^i^ijk I gijgik- 



(9) 



(10) 



Note that the tensor tjk is symmetric. Its trace is a sum 
of squares, and therefore always non-negative: 



tkk — gikgik > 0, 



(11) 



where we used the vanishing of partial traces, Eq. (|3J| . We 
note that this gives the well known result that U cannot 
have a local maximum in empty space, since a maximum 
would imply three negative eigenvalues of t and thus a 
negative trace. The absence of maxima in empty space is 
known as Wing's theorem ^^ and is here retrieved by a 
different route. Note that the non-negative trace relies on 
B being irrotational and divergence free. Therefore the 
same conclusion holds for electrostatic fields in vacuum. 
On the other hand time-dependent E and B fields are 
not irrotational and in fact do allow for a maximum of 
field magnitude in empty space |l6lll7|. 



C. Trapping frequencies 



f/: 



iUi + 2uigijXj + {uiCi-jk + gi]gtk)xjXk + h.o.t. (7) 



For a stationary point in Xi = we set dpU — 2uigip — 0. 
A trivial solution is that the field is zero, Ui = 0. Note 
however that this is generally only a stationary point of 



by: 



The potential for an atom in a magnetic field is given 



V = mpgFtJ'BVU, 



(12) 



with mp the magnetic quantum number, gp the Lande 
factor and jib the Bohr magneton. If we make a harmonic 
approximation around the potential minimum we find 
that the trap frequencies are given by: 



mpgrf^BTn 



(13) 



with u = ^UiUi, with T„ the eigenvalues of tjk and m the 
mass of the atom. Here we assume that the eigenvalues 

r„ > 0. 



Combining this expression with Eq. Hll|l we find 



u 



(14) 



Thus, remarkably, we find that this combination of trap 
frequencies is independent of the curvature Cijk and de- 
pends only on the gradient gij and the uniform field Ui. 



D. Where can a IP trap be created? 

Having established that stationary points, including 
IP traps, can only be found on the manifold A^, we now 
address the question whether a IP trap can be created in 
any arbitrary point on the manifold M . For convenience 
we choose a coordinate frame that diagonalizes gij , such 
that the zero eigenvector lies along coordinate direction 
63. The gradient tensor then takes a very simple form, 
with gii = —g22 = a as the only nonzero components. 
Since Ui must be chosen along the zero eigenvector, we 
write u = M3 63. We can then write the tensor tjk as 




+ '«3C3jfe- 



(15) 



Thus tjk depends on a single parameter M3 that is a mul- 
tiplier of the symmetric and traceless tensor c^jk- 

We can easily see the qualitative behavior of the eigen- 
values of tjk as a function of U3. Obviously, for U3 = 
the eigenvalues are (Ti,r2,T3) = (a'^,a^,0). For small 
values of u^ we can use perturbation theory to obtain 
the lowest eigenvalue to first order, yielding 



T. 



U3C333- 



(16) 



This shows that T3 can be made either positive or nega- 
tive by choosing the sign of u^. For small values of M3 the 
other two eigenvalues will remain positive. This means 
that for small U3 the stationary point will be either a IP 
trap or a saddle point with Morse index of 1, Morse index 
being the number of negative eigenvalues. Note that we 
can only make a IP trap if C333 y^ 0. In fact, for the case 
that C333 = a counterexample is easily found. 

For large enough positive or negative values of U3 the 
term c^jk will become the dominant term. Since c^jk is 



traceless, it has signature (-|-,-|-,— ) or (-1-,—,—). There- 
fore for large values of U3 we will always have a saddle 
point. The sign of U3 will determine whether the Morse 
index is 1 or 2. 

Since for small W3 the eigenvalues of tjk are given by 
(a^ -I- M3C1 13 , a^ -t- M3C223 , W3 C333 ) , we can tune the two non 
axial trap frequencies with u^ but the axial frequency is 
fixed by C333 [Eq. 13)]. 



IV. MOVING STATIONARY POINTS 

We now investigate how stationary points can be 
moved over the manifold by changing the uniform field 
Ui. We consider the situation that the spatial dependence 
of the magnetic field is defined, e.g., by a configuration of 
permanent magnets. We can influence the magnetic field 
pattern by applying a uniform external field. In terms 
of the above quantities, gij and Cijk are fixed, Ui is our 
control parameter. 



A. Moving loffe-Pritchard traps 

The movement of IP traps, which must clearly be con- 
strained to the manifold, is important in applications 
that require trapped atoms to move, such as beam split- 
ters and conveyor belts [l3|- We now calculate a dis- 
placement tensor djq = dxj/duq in a stationary point 
Xj = on the manifold. This tensor describes how the 
position of a stationary point moves when the uniform 
field is changed. To calculate it, we use the condition for 
a stationary point, dpU = 0, with U as in Eq. (TJ, 

dpU = 2uigip + 2{uiCijp -f gijgip)xj = 0. (17) 

Taking the derivative with respect to Uq and setting Xj = 
0, we solve for djq and find: 



_ dxj _^ 



(18) 



where {t~^)jp denotes the inverse tensor of tjp. In the 
basis where gpq is diagonal we see that djs — 0, i.e., a 
small field in the axial direction of a IP trap will not 
displace it. Since this means that the eigenvalue D^ 
is zero, djq is singular. Therefore djq is a mapping of 
a three-dimensional vector Uq onto a two-dimensional 
space, namely the tangential plane to the manifold. It 
is spanned by the two eigenvectors (Di,D2) correspond- 
ing to the nonzero eigenvalues. The vector Di x D2 is 
normal to the manifold and is found to be proportional 

to (ci33,C233,C333)- 

Note that during the movement the radial trap fre- 
quencies can be controlled using a bias field in the axial 
direction of the IP trap [Eq. (jT^ J. 



V. THE MANIFOLD AND FIELD ZEROES 

The manifold M is not only a powerful concept in the 
description of stationary points, it also has significance in 
the occurrence and movement of field zeroes. Such points 
of zero field are also minima of B and can thus serve as 
atom traps. However, these so-called quadrupole traps 
(QT) suffer from higher trap loss rates due to Majorana 
spin flips near the region of zero field. The movement of 
QTs is important in loading procedures, i.e. to transport 
atoms into the final IP trap. In this section we investigate 
how field zeroes can be moved and what happens when 
they approach the manifold. 

The manifold is the boundary between two regions of 
space y+, V^ , where det(gy) > and < respectively. 
Since g^ is traceless, and det(gij) is the product of the 
eigenvalues, gij must have two negative and one posi- 
tive eigenvalues in V~^ . Similarly, in V~ it has one neg- 
ative and two positive eigenvalues. It is impossible to 
move between the two regions V^ , V^ without one of 
the eigenvalues going through zero. This means that the 
manifold imposes restrictions on the movement of field 
minima that have a field zero. 



A. Moving quadrupole traps 

Since a quadrupole trap is a zero of the magnetic field, 
it is straightforward to give a prescription for how to 
move it by applying an external field. We call the sta- 
tionary field Bstat(r) and the desired trajectory r{t). In 
order to move the field zero along the trajectory, all we 
need to do is to use the external field Boxt(i) to cancel 
the local magnetic field. 



BextW = -B,tat(r(i)). 



(19) 



Like for the movement of stationary points, we can 
also express the movement of zeroes in terms of a dis- 
placement tensor dxj/duq, Eq. H18|l . For field zeroes this 
tensor takes a very simple form. 



= -i9-')M- 



(20) 



Thus, generally speaking, field zeroes do not disappear 
when the uniform field Uj is changed, they simply move 
through 3D space. 

The situation is different when a field zero approaches 
the manifold between the regions V^ , V^ . On the man- 
ifold det{gij) = so that djq does not exist. In fact, 
since one of the eigenvalues of gij vanishes, the magnetic 
confinement along one direction of the magnetic QT van- 
ishes. 



B. Pairs of zero field 

Having established how field zeroes (QTs) and IP traps 
move, the natural question arises what happens when a 



field zero approaches the manifold. We have just noted 
that field zeroes can be broadly categorized in two types, 
according to whether the signature of the gradient is 
(+,-!-,—) or (-1-,—,—). It turns out that the approach 
of the manifold by a field zero is accompanied by the ap- 
proach by another field zero, of the opposite type, from 
the other side of the manifold. 

If a zero is sufficiently close to A^ , we first assume that 
we can choose a point Xi = on Al such that the position 
of the zero Xi = £,i is in the direction G3 of the local 
quadrupole axis. Thus we have S.i — S.2 — and ^3 7^ 0. 
Furthermore g^ = on Al so that the requirement of 
zero field in S,i simplifies to 



C.33C3 = 0- 



(21) 



This shows immediately that the replacement ^3 — > —1^3 
yields another zero. The two zeroes are symmetrically 
placed around the point on At, in the direction of the 
G3 axis. The zeroes must be of opposite type, since 
they are on opposite sides of Ai. Finally, we note that 
the solution for ^3 of the above Eq. (|21|) . namely ^3 = 
zty^—Ui/ci33 implies that the local field direction m must 
be proportional to 0^33. This is just the normal vector to 
the manifold as mentioned above, so the local field Ui is 
normal to the manifold A4. 

The choice of a point Xi — on M such that ^i \\ G3 
is possible unless G3 is parallel to the manifold. In this 
special case, where C333 = 0, a zero arriving on the man- 
ifold can transform into a line of zero field, which lies 
entirely on the manifold. Such lines of zero require a 
treatment that goes beyond the second order field ex- 
pansion. Evidence from numerical examples and some 
highly symmetric analytical examples suggests that such 
lines are closed loops on the manifold. Small perturb- 
ing fields can split this loop into a number of zero pairs, 
which can be macroscopically separated. 

Finally we can also address the question whether the 
merging of field zeroes must necessarily take place on M . 
We assume that two field zeroes are sufficiently close so 
that it is sufficient to expand the field upto second order. 
We define local coordinates Xi around the first zero and 
x^ around the other. For the field expansion we can then 
write 



1 



JijXj 



ey^ijk-'^j-^k 



(22) 



Note that both Ui and u'^ are zero because we have two 
field zeros. We introduce the separation vector ^i and 
substitute Xi — S,i + x'^. Equating terms of equal powers 
in x'l we then obtain {gij -f g'ij)£,j = and thus Aet{gij + 
g[A — 0. Thus, in the second order field expansion, the 
midway point between the two zeroes Xj and x' lies on 
the manifold. Furthermore, the two field zeroes must be 
of opposite type. Apparently, in order for field zeroes of 
similar type to approach each other the leading term in 
the field expansion must be higher than second order. 




FIG. 1: (Color online) Manifold corresponding to the stan- 
dard loffe-Pritchard trap, described by the field of Eq. 12311 
using e = 1. The black bars indicate the orientation of the 
loffe bars. 




VI. PRACTICAL lOFFE TRAPS 

In this section we consider the shape of the manifold 
det(gij) = in some cases of practical interest. 

A. Standard loffe-Pritchard trap 



FIG. 2: (Color online) Manifold created by a Z-shaped wire. 
Since the manifold is given by the gradient its shape does not 
depend on the bias field. The current is only a multiplier 
and is also of no importance for the shape. The vectors Di, 
D2, G3 and Di x D2 have been drawn in the point (0,0,1). 
The vector G3 gives the axial direction of the IP, Di x D2 is 
normal to the manifold. 



The prototypical IP trap [ifl UM consists of four long 
current carrying wires ( "lofFe bars" ) , for creating a cylin- 
drical quadrupole field, in combination with two pinch 
coils creating confinement in the axial direction. The 
field at the IP can be approximated by: 



(23) 



Where u, a and c are the uniform field, radial gradient 
and axial curvature, respectively. The manifold produced 
by this field is shown in Fig. ^ It is described by the 
equation: 




e(x2 - y'^) + z{2z'^ + y^ ■ 



2^] 







(24) 



Where e = 2a/ c. Note that for e = this describes 
the flat surface 2 == 0. The hole in the manifold has 
typical size e. In the point (0, 0, 0) the eigenvectors of djq 
point in the cc, y, and z directions and so the lowest order 
displacement of the IP under influence of an external field 
is in a direction prependicular to the axial direction. The 
shape shown is only realistic in the region where Eq. H23|l 
is a good approximation of the field. 




FIG. 3; (Color online) The azimuthal angle tp of G3 with 
respect to the x-axis as a function of 6. Note that directly 
above the middle of the central wire G3 always lies in the 
xy-p\a.ne. The inset shows the geometrical meaning of 6. 



B. Z-wire loffe trap 

A method routinely used in atom chip experiments for 
creating IP traps involves a current carrying wire bent in 
a Z-shape in combination with a uniform bias field [l9j. 
In Fig. 12 the manifold for such a wire is shown together 
with the eigenvectors of djq and the IP axis G3. The 
vector G3 straight above the middle of the central wire 



always lies in the xy-plane. We find that the azimuthal 
angle ip between G3 and the x-axis is given by: 



<P 



arctan 



cos 26 sec 9 
2 + cot^ 6 



(25) 



Here tan6' = 2z/s, where s is the separation between 
the two wires in the y-direction and z is the height above 
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FIG. 5: (Color online) Two-dimensional cross-section showing 
the magnetic potential and the manifold at a = 10/im. The 
IP traps can be transferred over the manifold. For one IP the 
route to the neighbouring lattice site is shown. The arrows 
indicate how the trap axis varies along the way. 



FIG. 4: (Color online) Array of magnetic material (dark re- 
gions). The magnetization is out of plane. The equivalent 
edge current forms "z" -shapes at every lattice point. In com- 
bination with a bias field in the y-direction IP traps can be 
created. The shape of the manifold is shown in the inset. 



the central wire. This result is plotted in Fig. |31 For 
z — 0, Gs points in the a;-direction, in the limit z -^ oo 
it points in the y-direction. For a given s, G3 points in 
the x-direction if we choose z ~ ^s- 



C. Array of lofFe traps 

The exact shape of the manifold is of particular im- 
portance in arrays of IP traps that might be used as 
conveyor belts/shift registers. Such devices are promis- 
ing for quantum information processing applications |2Q| . 
Atoms sitting at a lattice site that is connected to another 
via the manifold can be shifted there using a uniform bias 
field, while remaining a IP trap. Moreover, the trap fre- 
quencies along the way can be tuned using Eq. (|13|l . 

To make these ideas more explicit, we discuss an array 
of IP traps created by permanent magnetic material in 
combination with a uniform bias field. In Fig.^the array 
of magnetic material is shown. Since its magnetization 
is out of plane we can think of the material as having an 
equivalent current of magnitude Mh running around its 
edges, where M is the magnetization and h is the height 
of the material. This equivalent current forms a Z-shape 
at every lattice site in the array. We apply a uniform bias 
field in the y-direction to create loffe traps at all lattice 



sites. 

In this particular design the magnetization is 
M =800 kA/m and the height of the material is 
h =250 nm. A bias field of 17.9 G in the y-direction 
then produces IP traps at 10 /im from the surface with 
trap frequencies (21, 20, 5.4) kHz and a residual field of 
4.5 G at the trap bottom. 

We are interested in where the individual IP traps can 
be moved. Therefore we draw the manifold over a region 
of the array containing several unit cells. As can be seen 
in Fig. ^ IP traps are only connected in a diagonal di- 
rection from lower right to upper left. This means that 
this array can be used as a shift register only in this di- 
rection. The manifold clearly does not allow shifting the 
IP traps in the perpendicular direction. One could try to 
move the traps in the perpendicular direction as field ze- 
roes (QT). However, we find that this leads to a sequence 
of splitting and recombination of zeroes, every time the 
manifold is crossed. Thus, in spite of its appearance, this 
is a ID shift register. 

Finally, for the sake of completeness, we show in Fig. [51 
how the axial direction of the loffe traps varies over the 
manifold. 



VII. CONCLUSION 

In conclusion, we have shown that loffe Pritchard traps 
as well as other stationary points of B(r) are confined 
to a two-dimensional curved manifold M defined by 
det{gij) = det{dBi/dxj) = 0. Furthermore, in any point 
of Ai where the local quadrupole axis is not parallel to 
A^ a IP trap or other stationary point can be created by 



choosing the magnetic field paraUel to this axis, i.e. the 
eigenvector of gij corresponding to eigenvalue zero. We 
have given an expression for the movement of stationary 
points over the manifold, in response to a change of an 
external uniform control field. 

We have shown that the same manifold plays an im- 
portant role in the behavior of field zeroes, and separates 
regions of space where field zeroes are of an opposite type. 
These field zeros can be moved anywhere in their respec- 
tive region of space but can only disappear on the mani- 
fold. In this annihilation two field zeros of opposite type 
merge and form a loffe Pritchard trap. Similarly, a IP on 
the manifold can be made to split into two field zeroes of 
opposite type, on opposite sides of the manifold. 

The manifold is a new conceptual tool in designing 
magnetostatic or electrostatic potentials. Understanding 
how IP traps can be moved over this manifold is of great 
use in experiments that require some level of dynamics 
such as shift registers for trapped atoms or double well 



experiments. The shape of the manifold is important in 
loading protocols for permanent magnetic atom chips, 
defining possible routes for transfering loffe Pritchard 
traps. 
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